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Many-body localization characterized from a one-particle perspective 
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We show that the one-particle density matrix p can be used to characterize the interaction-driven 
many-body localization transition in closed fermionic systems. The natural orbitals (the eigenstates 
of p) are localized in the many-body localized phase and spread out when one enters the delocalized 
phase, while the occupation spectrum (the set of eigenvalues of p) reveals the distinctive Fock- 
space structure of the many-body eigenstates, exhibiting a step-like discontinuity in the localized 
phase. The associated one-particle occupation entropy is small in the localized phase and large in the 
delocalized phase, with diverging fluctuations at the transition. We analyze the inverse participation 
ratio of the natural orbitals and find that it is independent of system size in the localized phase. 

PACS numbers: 72.15.Rn 05.30.Rt 05.30.Fk 


Introduction. While the theory of noninteracting dis¬ 
ordered systems is well developed mm, the possibility of 
a localization transition in closed interacting systems has 
only recently been firmly established I5T-E5], This many- 
body localization (MBL) transition occurs at finite en¬ 
ergy densities and is not a conventional thermodynamic 
transition [23] 125] . Instead, it can be understood as a 
dynamical phase transition, associated with the emer¬ 
gence of a complete set of local conserved quantities in 
the localized phase, which thus behaves as an integrable 
system (261(301 . This restricts the entanglement entropy 
of the eigenstates to an area law [3T] , in contrast to the 
volume law predicted by the eigenstate thermalization 
hypothesis for the ergodic delocalized phase [32H53] . At 
the localization transition, the fluctuations of the entan¬ 
glement entropy diverge [[16] [35]. The effects of MBL 
are also observed in the dynamics following, for example, 
a global quench from a product state, wherein dephas¬ 
ing between the effective degrees of freedom leads to a 
characteristic logarithmic growth of the entanglement en¬ 
tropy These features comprise a much richer 

set of signatures than in the context of noninteracting 
systems, for which, in the spirit of one-parameter scaling, 
the notion of a localization length based on single-particle 
wave functions generally suffices mm- 

In view of the rich phenomenology of many-body local¬ 
ization it is natural to ask, both from a fundamental point 
of view as well as for the interpretation of experimental 
data |36j 37], to which extent (if at all) the MBL tran¬ 
sition can be detected and characterized from a single¬ 
particle perspective. Here we show that a rather com¬ 
plete characterization of many-body localization is in¬ 
deed possible based on the eigenvalues (occupations) and 
eigenstates (natural orbitals) of the one-particle density 
matrix. The one-particle density matrix was originally 
introduced by Onsager and Penrose to extend the notion 
of a Bose-Einstein condensate to interacting systems [3S]. 


Importantly, the natural orbitals take a Bloch form in 
translationally invariant systems, providing a true many- 
body generalization of the Bloch theorem [35]. This nat¬ 
urally suggests studying the effect of disorder, as in recent 
studies of localization of thermalized hard-core bosons in 
a quasi-periodic potential !40, EH]. However, so far no 
connection to many-body localization has been made. 

We are further motivated to consider the one-particle 
density matrix because it naturally focusses on the dy¬ 
namics of one particle in the presence of all the others, 
without the need to resort to a mean-field theory or to 
sacrifice particle indistinguishability. As we will see, this 
perspective retains sufficient information to capture the 
genuine many-body aspects that set many-body local¬ 
ization apart from Anderson localization transitions in 
noninteracting systems. In particular, the occupations 
sharply reorganize themselves from being close to either 
zero or one in the localized phase to being in between 
these extremal values in the delocalized phase, thus re¬ 
flecting a delocalization transition in Fock space that cor¬ 
responds to a mixing of product states. It follows that 
in the localized phase the occupation spectrum devel¬ 
ops a step-like discontinuity, similar to a Fermi-liquid. 
The associated one-particle occupation entropy is large 
and proportional to the system size in the delocalized 
phase, corresponding to the volume law of thermal states, 
while in the localized phase it is small. The dynamics of 
one particle in the effective bath of the others thus pro¬ 
vides complementary information to the dynamics of a 
spatially confined region in the effective bath of its sur¬ 
rounding. In addition, we show that the transition leaves 
direct signatures in the natural orbitals, which are local¬ 
ized in the many-body localized phase and spread out 
over the system when one enters the delocalized phase. 
We show that the inverse participation ratio (IPR) of the 
natural orbitals depends on the system size in the delo¬ 
calized phase, while it is independent of system size in 
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the localized phase. 

Model and method. We consider spinless fermions in 
one dimension with a nearest-neighbor repulsion and di¬ 
agonal disorder, described by the Hamiltonian 



Here c\ creates a fermion on site i = 1,2,... ,L and 
Hi = c\ci is the associated number operator. Energies 
are expressed in units of the nearest-neighbor hopping 
constant t, so that V is a dimensionless measure of the 
strength of the nearest-neighbor repulsive interactions. 
The diagonal disorder is introduced via a box distribu¬ 
tion of the onsite potentials £ [—IT, W], We study 
this system using exact diagonalization at finite sizes 
L = 10,12,14 (10 5 disorder realizations), L = 16 (10 3 
realizations) and L = 18 (500 realizations), imposing pe¬ 
riodic boundary conditions and fixing the overall occu¬ 
pation at half filling (number of particles N = L/ 2). We 
mainly focus on the energy region around the band cen¬ 
ter, e = 1 where e = 2 (E - E min ) / (E max - E min ) with 
Bmax and E m - m the maximum and minimum energy for 
each disorder realization, and take the 6 eigenstates clos¬ 
est to this energy. This energy corresponds to infinite 
temperature in the thermodynamic limit. At the fixed 
interaction strength V = 1, the critical disorder strength 
W c is found to take values in the range between three 
and four mmmm- 

Given a many-body eigenstate \ip n ) of the Hamilto¬ 
nian 0. the one-particle density matrix is defined as 

Pij = W’nlcJcjIV’n)- (2) 

The natural orbitals \<j> a ) with a = 1,2,... ,L, are ob¬ 
tained by diagonalizing p , 

p\(f> a ) = n a \(f> a ), (3) 

which delivers a basis of single-particle states. The 
eigenvalues n a are interpreted as occupations, with 
'}2a=i n a = tr p = N equal to the total number of par¬ 
ticles in the system. We order the natural orbitals by 
descending occupation, n\ > n ,2 > ... > ul- 

Occupation spectrum. In a noninteracting fermionic 
system, barring degeneracies, each many-body eigenstate 
\ipn) can be written as a Slater determinant of N single¬ 
particle states. The occupations from the one-particle 
density matrix are then fixed to n a = 0 or 1, with the 
natural orbitals with n a = 1 spanning the space of the 
single-particle states used in the Slater determinant. Oc¬ 
cupations departing from n a = 0 or 1 can therefore be 
interpreted as a signature of the true many-body nature 
of the eigenstates in the interacting system and, thus, 





FIG. 1. (Color online) (a) The main panel shows the disorder- 
averaged occupation of the natural orbitals with index a for 
different values of disorder strength (system size L = 16, in¬ 
teraction strength V = 1). The dotted line shows the oc¬ 
cupation in a noninteracting system, which is independent 
of the disorder strength. The vertical line indicates the ex¬ 
pected discontinuous behavior of the occupations in the in¬ 
finite system-size limit of the MBL phase. The emergence 
of this discontinuity is further illustrated in the inset, which 
shows the average occupations (njv) and (njv+i) as a function 
of disorder strength, for L = 16 and three values of the inter¬ 
action (V = 0.5,1,1.5). (b) Occupation discontinuity An as 
a function of disorder strength and energy density (left panel) 
or interaction strength (right panel), (c) Distribution of oc¬ 
cupations n a in the delocalized phase (W = 0.4,0.8), near 
the MBL transition (W = 3.0), and deep in the localized 
phase (W = 6.0,8.0), for system size L = 16 and interaction 
strength V — 1. 


as a proxy of the delocalization of such states in Fock 
space. We expect this Fock-space delocalization to be 
pronounced in the metallic phase, while it should be sup¬ 
pressed in the MBL phase M. 

In Fig.JIJa) we show the disorder-averaged occupations 
for different values of disorder for L = 16, with (.) denot- 
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ing the disorder average. The horizontal dashed lines rep¬ 
resent the occupations (n a ) = 1 for a < N, (n a ) = 0 for 
a > N in the noninteracting limit V = 0, where the sys¬ 
tem is Anderson localized for any finite disorder strength. 
The quasi-discontinuous jump An = njv — njv+i = 1 
between these values is indicated by a vertical line. 
In the interacting system, deep in the localized phase 
(W = 6,8), half of the natural orbitals remain almost 
fully occupied, (n a ) ss 1, with the other half being al¬ 
most unoccupied, (n a ) ~ 0. As one decreases the disor¬ 
der and approaches the transition {W = 3), more orbitals 
acquire a finite occupation, while for even smaller disor¬ 
der, in the delocalized phase (W = 0.4,0.8), the occupa¬ 
tion of all orbitals approaches the mean filling fraction, 
(n a ) » N/L = 1 / 2 . 

The redistribution of occupations with decreasing dis¬ 
order goes along with a reduction of the step-like be¬ 
havior quantified by An. A more detailed view of this 
aspect is provided by the inset of Fig. [lja) , showing the 
disorder dependence of the occupations (tin), (jijv+i) for 
three values of interaction strength V = 0.5,1.0,1.5. In 
the delocalized phase, both occupations are close to the 
mean filling fraction, (n a ) ~ N/L = 1/2, while deep in 
the localized phase they tend to their asymptotic val¬ 
ues (tin) = 1, (n/v+i) = 0 45]. The dependence of 
the discontinuity An on energy density, shown in the 
left panel of Fig. TVb) , recovers the many-body mobil¬ 
ity edge 0 QH mrsl], while at small and large interac¬ 
tion strengths (right panel) the delocalized phase shrinks, 
consistent with observations from dynamics in the same 
model [?3]. According to these results, the occupation 
spectrum serves as a reliable indicator of many-body lo¬ 
calization. 

One-particle occupation entropy. A well documented 
aspect of MBL is the appearance of strong fluctuations 
around the localization-delocalization transition uadi 
SM3- In terms of the occupations, this is addressed in 
Fig. s c), which shows the probability distribution func¬ 
tions P(n a ) for different disorder strengths in a semi¬ 
log plot. In the large disorder limit the distribution 
is bimodal with peaks at n a = 0,1, with very little 
weight in the central region between these extremal val¬ 
ues. This bimodal distribution is characteristic of the lo¬ 
calized state, in analogy to the noninteracting scenario. 
As expected, close to the transition (W = 3) the distri¬ 
bution is wide, with significant weight across the whole 
range of occupations. Finally, in the delocalized phase 
with low enough disorder the distribution becomes again 
narrower, but now is concentrated around the filling frac¬ 
tion N/L = 1/2. 

In order to quantify these fluctuations further, we con¬ 
sider the entropy 

S = —trplnp = — ^ n a ln(n Q ). (4) 

OL 

As this entropy is determined by the occupations of the 



FIG. 2. (Color online) Dependence of the disorder averaged 
one-particle occupation entropy (S), defined in Eq. , on 
the disorder strength, for different system sizes at interaction 
strength V = 1. The dashed lines indicate the maximal value 
L In 2, corresponding to the volume law for the entropy in a 
fully delocalized system. In contrast, in the localized phase 
the entropy becomes small. Inset: Variance var S of the en¬ 
tropy due to sample-to-sample fluctuations in the disorder 
ensemble as a function of disorder strength, for different sys¬ 
tem sizes at interaction strength V = 1. The peak in the 
variance indicates the location of the MBL transition. 


natural orbitals we call this the one-particle occupation 
entropy, to distinguish it from the entanglement entropy 
of the many-body eigenstates. The disorder-averaged en¬ 
tropy (S) is shown in the main panel of Fig. [2j as a func¬ 
tion of disorder strength for different system sizes. In 
the delocalized phase the entropy approaches the max¬ 
imal value lln2, indicated by the dashed lines. This 
corresponds to a volume law as displayed, in general, 
by extensive thermodynamic properties and many-body 
eigenstates in ergodic systems. In contrast, the entropy 
in the localized phase is much smaller (35]. 

In the inset of Fig. [2] we show the variance var S = 
(S 2 ) — ( S ) 2 of the entropy as a function of disorder 
strength. For the corresponding case of the entanglement 
entropy, it is known hue] that the variance vanishes 
in the thermodynamic limit (L —> oo) both in the lo¬ 
calized and in the delocalized phase, where in the latter 
phase this is consistent with the eigenstate thermaliza- 
tion hypothesis. Furthermore, in finite systems, the vari¬ 
ance of the entanglement entropy is sharply peaked in 
the crossover regime, which is associated with the mixing 
and coexistence of localized and delocalized regions near 
the transition, becoming sharper with increasing system 
size [16] . This universal behavior of the entanglement en¬ 
tropy is mirrored by the one-particle occupation entropy. 
The occupation spectrum therefore recovers a reliable sig¬ 
nature of the MBL transition, giving quantitative access 
to the locus of the transition. 

Delocalization of natural orbitals. Because of the 
emerging degeneracy of the occupations deep in the local¬ 
ized and delocalized phase, one may wonder whether the 
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FIG. 3. (Color online) Evolution of the probability distribu¬ 
tion of the IPR for increasing system size (a) in the delocalized 
phase (W = 0.4), (b) close to the transition {W = 3.0), and 
(c) deep in the localized phase (W = 8.0). The insets in 
(a) and (c) show examples of the natural orbitals, (d) Av¬ 
erage participation ratio (1/IPR) = £ of the natural orbitals 
as a function of disorder strength. In the localized phase £ 
is independent of the system size, while for small disorder 
it saturates at £ ss L/2. Inset: Average participation ratio 
as a function of interaction strength V for several values of 
disorder strength (L = 16). 

natural orbitals themselves display any signatures of the 
MBL transition. As we show in the insets of Figs. [3j)a) 
and (c), the orbitals indeed turn out to be well localized 
in the MBL phase, while they are far more extended in 
the delocalized phase. The multiply-peaked structure of 
the natural orbitals for weak disorder suggests that the 
delocalization transition involves the formation of a chain 
throughout the system via which the particle can hop 
resonantly. Given the complete set of natural orbitals, 
a measure of the localization of the occupied states can 
then be derived from the inverse participation ratio 

IPR= ^E ri “Ei^wi 4 - (5) 

a =1 i —1 

The IPR is normalized to take the maximal value 1 for 
a system in which all occupied states are fully localized, 
while it takes the minimal value 1/L when all occupied 
states are fully extended. In between these two extremes, 


the resonant-hopping picture for the delocalized phase 
suggests that the IPR scales inversely with the system 
size, while in the localized phase it should be indepen¬ 
dent of system size. These tendencies are confirmed in 
the main panels of Figs. [3ja-c) , which show, for three 
disorder strengths in the delocalized, transitional, and 
localized regime, how the probability distribution of the 
IPR depends on the system size. In the delocalized phase 
(a), the flow with system size is indicative of a 1/L be¬ 
havior, while in the localized phase (c) the distribution 
is almost independent of system size, with a peak close 
to the maximal value IPR = 1. Close to the transi¬ 
tion (b), the IPR distribution is wide, with no discernible 
trend with system size. It is therefore suggestive to intro¬ 
duce the characteristic length £ = (1/IPR). Figure [3](d) 
shows the disorder-strength-dependence of £ for different 
system sizes. In the localized regime this characteristic 
length is independent of system size. With decreasing 
disorder strength £ increases, whereas at very small dis¬ 
order it approaches the value £ ~ L/2. While £ is still 
small at the transition in the accessible system sizes, the 
orbitals spread out significantly once one enters into the 
delocalized phase. Moreover, as shown in the inset of 
Fig- §d), £ depends non-monotonically on V: it first 
increases as V increases, then takes a maximum at a In¬ 
dependent value and finally decreases again in the large 
V-limit. A similar behavior was observed in spectral fluc¬ 
tuations in this model [33] . Our quantity £ thus captures 
the delocalizing effect of interactions, both in the delo¬ 
calized and in the MBL phase. 

Summary and outlook. In conclusion, the one-particle 
density matrix uncovers essential many-body aspects of 
interacting disordered fermions. Our results suggest that 
in the thermodynamic limit the one-particle occupation 
spectrum is continuous in the delocalized phase but de¬ 
velops a finite discontinuity in the localized phase. The 
corresponding occupation entropy shares features with 
the many-body entanglement entropy, one of the prin¬ 
cipal vehicles for the theoretical characterization of the 
many-body localization transition. The delocalization is 
also observed in the structure of the natural orbitals, 
which is reflected in a system-size dependent inverse par¬ 
ticipation ratio. These findings support the conceptual 
picture that the many-body localization transition in¬ 
volves delocalization both in Fock space and in real space. 
Our approach should therefore apply to a broad range of 
systems that follow this scenario, which can be further 
enriched when the particle number is not conserved. An 
interesting and timely application of our work would con¬ 
sist in analyzing the one-particle density matrix for the 
system that was experimentally realized in Ref. 1371 
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Appendix A: System-size dependence of the 
occupation spectrum 


In the main text, we showed that the occupation spec¬ 
trum ( n a ) derived from the one-particle density matrix 
allows to distinguish the many-body localized phase from 
the delocalized phase. In the localized phase, the occu¬ 
pations tend to be close to 1 and 0, which suggests a 
discontinuous dependence (finite jump An) to persist in 
the thermodynamic limit. Likewise, in the delocalized 
phase, one expects that the eigenstate thermalization hy¬ 
pothesis holds, which corresponds to a smooth behavior 
of the occupation spectrum in the thermodynamic limit. 
While the system sizes accessible in exact numerics are 
still relatively small, Fig. [I] shows that these expecta¬ 
tions are consistent with our numerical results. In the 
delocalized phase [panels (a)-(c)], increasing the system 
size smoothes out the occupation spectrum (closer to the 
transition this happens more slowly), while in the local¬ 
ized phase [panels (d)-(f)] we maintain two branches at 
occupations close to 1 or 0. In a more detailed scenario 
still consistent with this data, the jump An would con¬ 
tinuously increase from 0 at the MBL transition to 1 deep 
in the MBL phase. 
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FIG. 4. Occupation spectrum ( n a } at fixed disorder strengths for (a) W = 0.4, (b) W = 0.8, (c) W = 2.0, (d) W = 3.0, (e) 
W = 6.0, and (f) W = 8.0, for different system sizes L [cf. main panel of Fig. 1(a) in the main text]. The index a is scaled so 
that in the thermodynamic limit the horizontal axis runs from 0 to 1. The interaction strength is V = 1. 



FIG. 5. System-size dependence of the occupation entropy (S) 
versus L for disorder strength W = 0.4, 0.8, 2.0,3.0,6.0, 8.0 
with interaction strength V = 1. 


the overall trend discussed in the main text: S is large in 
the delocalized phase and then becomes much smaller in 
the many-body localized phase, approaching zero for very 
strong disorder. In the delocalized phase, the entropy is 
proportional to L and becomes independent of disorder 
at large system sizes, consistent with thermalization. For 
the accessible system sizes, we observe that the entropy 
still increases in the localized phase, which is related to 
the existence of occupations n a ^ 0,1. In fact, the en¬ 
tanglement entropy also has a residual L-dependence for 
values of W in the localized phase close to the transi¬ 
tion m- As a caution, we stress that the accessible 
system sizes may be too small to resolve the asymptotic 
A-dependencies of either entropy for values of W suffi¬ 
ciently close to the transition. As long as the localization 
length is comparable to the system size, we do not expect 
to be able to resolve an area law numerically. 


Appendix B: System-size dependence of the 
occupation entropy 


Figure [5] shows the L-dependence of the occupation en¬ 
tropy for several values of L. The data is consistent with 










